x1=[1.4	2	2	3.3	7	6	4.7	3.5	2.5	2.3	3	3	2	3	3	1.9	2	5	2	1.7	1.3	1	1.7	1	2.3	3	2	4	4	2.8	4	3.3	2.9	3.2	4	3.7	3	2.9	3	4.1	3	4	2	1.9	2	1.9	2	3	3]';
x2=[1.4115	1.9913	2.6401	3.1966	6.6236	7.4063	4.9308	3.8367	2.7346	1.9593	1.927	2.8445	2.5082	1.6078	2.8519	2.5356	1.4744	1.8695	4.2245	1.8107	1.4898	1.1373	0.8744	1.4338	0.9255	1.9834	2.5829	1.8785	3.4231	3.5202	2.7233	3.4791	3.0505	2.6982	2.915	3.5696	3.4195	2.851	2.6731	2.787	3.5796	3.0042	3.5293	2.2189	1.8627	1.8563	1.7848	1.8528	2.6114]';
x3=[1.3953	2.0204	2.3262	4.0457	5.1075	6.0042	5.4131	3.9433	2.989	2.2209	1.5686	1.7177	2.2169	2.2787	1.378	2.3908	2.1954	1.2866	1.5285	3.7826	1.5716	1.2962	0.9884	0.7597	1.2652	0.8075	1.7165	2.3339	1.6929	3.022	3.3661	2.4785	3.1844	2.8289	2.5603	2.7248	3.306	3.1776	2.624	2.5788	2.5678	3.4776	2.8066	3.3383	2.0972	1.7564	1.755	1.7045	1.7647]';
%x4=[0.3726	0.545	0.7812	1.0975	1.564	1.6607	2.5494	5.2061	5.0886	4.1956	3.2086	2.3357	2.0211	2.4237	2.5071	1.8562	2.3895	2.5	1.7879	1.7071	3.695	2.1025	1.5719	1.1946	0.9159	1.3259	0.9431	1.7324	2.3638	1.8265	3.0516	3.3055	2.5563	3.2029	2.8691	2.5327	2.6634	3.2251	3.141	2.6558	2.4884	2.5206	3.2623	2.681	3.2287	2.0058	1.6855	1.6859	1.6191]';
x5=[0.1761	0.3536	0.5173	0.7415	1.0416	1.4844	1.5762	2.4197	4.9412	4.8297	3.9821	3.0453	2.2169	1.9183	2.3004	2.3796	1.7618	2.2679	2.3728	1.6969	1.6202	3.507	1.9955	1.4919	1.1338	0.8693	1.2584	0.8951	1.6443	2.2435	1.7336	2.8963	3.1373	2.4262	3.0399	2.7231	2.4038	2.5278	3.0609	2.9812	2.5207	2.3618	2.3923	3.0963	2.5446	3.0644	1.9037	1.5997	1.6001]';
num=length(x1);
res=0
for i=1:20
   res=res+(abs(x5(i,1)-x1(i,1))/x1(i,1));
end
res/40
plot(1:num,x1,'r+-',1:num,x2,'b+:',1:num,x3,'bo:',1:num,x5,'bx:');